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Abstract. We inquire the statistical properties of the pair formed by the Navier-Stokes 
equation for an incompressible velocity field and the advection-diffusion equation for a scalar 
field transported in the same flow in two dimensions (2d). The system is in a regime of 
fully developed turbulence stirred by forcing fields with Gaussian statistics, white-noise in 
time and self-similar in space. In this setting and if the stirring is concentrated at small 
spatial scales as if due to thermal fluctuations, it is possible to carry out a first-principle 
ultra-violet renormalization group analysis of the scaling behavior of the model. Kraichnan's 
phenomenological theory of two dimensional turbulence upholds the existence of an inertial 
range characterized by inverse energy transfer at scales larger than the stirring one. For 
our model Kraichnan's theory, however, implies scaling predictions radically discordant 
from the renormalization group results. We perform accurate numerical experiments to 
assess the actual statistical properties of 2d-turbulence with power-law stirring. Our results 
clearly indicate that an adapted version of Kraichnan's theory is consistent with the observed 
phenomenology. We also provide some theoretical scenarios to account for the discrepancy 
between renormalization group analysis and the observed phenomenology. 



1. Introduction 



Two-dimensional (2d) turbulence is interesting for several reasons. In laboratory experiments 
2d turbulence has been realized and studied with electromagnetically driven liquid metals 
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||62l [6TI [T5l and thin soap films [|63l [641 [651 . In geophysical flows, vertical confinement 
suggests the possibility to describe the mesoscale dynamics of atmosphere and oceans in terms 
of two-dimensional fluid models 11201 |24"1 . Indeed, observational data such as the Nastrom- 
Gage spectrum 11571 [581 . studies based on the MOZAIC database ll55l and on the EOLE 
Lagrangian balloons in the low stratosphere |[43l support the existence of a mesoscale —5/3 
power-law energy spectrum which may be the consequence of a two-dimensional inverse 
cascade. Although recent studies Il46l |49l suggest the occurrence of a fairly more complex 
physical phenomenology (see e.g. Il69l 1701 ). the 2c/-approximation remains an important 
benchmark for understanding the atmospheric physics at synoptics and planetary scales P71 
as well as in other geophysical contexts (see ll671 l4T|| and references therein). For example, 
analysis of spectral kinetic energy fluxes in satellite altimeter data provides strong evidences 
of the occurrence of an inverse energy cascade in the ocean 14J. From the point of view of 
statistical mechanics, 2<i-turbulence is a prototype of non-equilibrium systems whose steady 
state is not described by Boltzmann statistics. 

At variance with its three-dimensional counterpart for which the total kinetic energy 
is the unique inviscid invariant, the two-dimensional Navier-Stokes equation preserves also 
the total enstrophy (in the absence of forcing and dissipation) 11301 . Enstrophy conservation 
is a key ingredient for the proof 11421 of the existence and uniqueness of the solution of 
the Cauchy problem for the 2d Navier-Stokes equation with deterministic forcing. Very 
recently, this result has been extended to stochastic stirring. In particular it was shown 
11401 [T6l l25l 1T71 1T81 l52l that the solution is a Markov process exponentially mixing in time 
and ergodic with a unique invariant (steady state) measure even when the forcing acts only on 
two Fourier modes [32J . 

A phenomenological theory proposed by Kraichnan [1371 and further extended by 
Batchelor and Leith 0H (see also [HH [6l d |46l ED) predicts the presence of a double 
cascade mechanism governing the transfer of energy and enstrophy in the limit of infinite 
inertial range. Accordingly, an inverse energy cascade with spectrum characterized by a 
scaling exponent —5/3 appears for values of the wave-number p smaller than pf, the typical 
forcing wave-number. For wave-numbers larger than pf a direct enstrophy cascade should 
occur. In this regime the energy spectrum should have a power-law exponent equal to —3 plus 
possible logarithmic corrections hypothesized in [|38l to ensure constancy of the enstrophy 
transfer rate. Very strong laboratory experiments reviewed in [[671 l36l . and numerical 
experiments (see e.g. |[P3l and references therein) corroborate Kraichnan's theory. A long 
standing hypothesis |[60l . which has recently found support through numerical experiments 
[El CGI, also surmises the existence of a conformal invariance underlying the inverse energy 
cascade of 2d turbulence. 

Despite these successes, a first principle derivation of the statistical properties of 2c? 
turbulence is still missing. An attempt in this direction has recently been undertaken 
lf34l [33l [35l El by inquiring the scaling properties of the velocity field and of the transported 
scalar field (passive scalar) when they are sustained by a random Gaussian forcing with self- 
similar spatial statistics. The Holder exponent e of the forcing correlation provides an order 
parameter interpolating between small-scale thermal stirring and large-scale stirring. 
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In three (and more) dimensions ultra-violet renormalization group analysis [|27l |23l of 
this model yields the result 1 — 4£:/3 for the scaling exponent of the kinetic energy spectrum 
holding to all orders in a perturbative expansion in powers of e. Kolmogorov scaling 11301 is 
recovered when the energy input becomes dominated by its infra-red components at e equal 
two. The results of [|27ll23l are coherent with physical intuition because only the case e equal 
two is a model of fully developed turbulence. Recent numerical simulations validate within 
available resolution such picture f[66l ITTTl . 

The extension of renormalization group analysis to the 2d case is instead not 
straightforward and was only achieved in [|34l . Although the prediction for the kinetic energy 
scaling exponent is the same as in the three dimensional case, the result cannot be easily 
reconciled with the phenomenological intuition based on Kraichnan's theory. The latter 
suggests the onset of an inverse energy cascade already at e equal zero when the energy and 
enstrophy input is dominated by the ultra-violet degrees of freedom. 

The purpose of the present work is to shed light on the apparent contradiction between the 
phenomenological and the renormalization group theory. In doing so, we extend and complete 
results presented in a previous letter 11531 . In particular, in section |2] we illustrate the details 
and interpretation of the model. In section [3] we summarize the results of the renormalization 
group analysis found in Il3"4ll33ll35l l3ll. In section @] we draw on [[6l [7]] to solve the Karman- 
Howarth-Monin equation under the hypotheses of Kraichnan's theory. In doing so, we focus 
first our attention on the parametric range where a direct comparison with the renormalization 
group theory is possible. We also shortly discuss the predictions of the phenomenological 
theory for the general case, as a benchmark for our numerical experiments. In section \5\ 
we compare the renormalization group predictions with the outcomes of a direct numerical 
integration of the model equations. There we give clear evidences that in all parametric ranges 
only the phenomenological theory a la Kraichnan is able to describe the observed behavior 
of the stochastic flow. Finally, in section [6] we discuss possible mechanisms underlying the 
discrepancy between the renormalization group predictions and the observed scaling behavior 
of the model. 

2. The model 

We consider the 2d Navier-Stokes equation governing the evolution of the velocity v of an 
incompressible Newtonian fluid 



where v is the kinematic viscosity and k the diffusivity of the scalar field. The Ekman friction 
term — - included in Eq.® ensures that a steady-state is attained by damping kinetic energy 
transfer towards larger and larger scales ll26ll . Both equations (OQ) and © are sustained by 



(d t + vd x )v = vd 2 x v-d x P-- + f 

T 

d x ■ v = d x ■ f = 



(1) 

(2) 



and the forced advection-diffusion equation for a scalar field 6. 
(d t + vd x )6 = Kd 2 x 6 + g 



(3) 
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stochastic forcing fields, respectively / and g, with Gaussian statistics such that 

-<f(x,t) y=^g(x,i) ^=0 (4) 

and 

-< f a (x, t)f(y, s) y= S(t - s) F al3 (x -y,m f ,M f ), a , (3 = 1, 2 (5) 

-< g(x, t)g(y, s) y= 5(t - s) G(x - y, m g , M g ) (6) 

Whilst time decorrelation in ©, © is meant to preserve Galilean invariance of the statistics 
(in the absence of Ekman friction), the spatial part of the forcing correlations is chosen to be 
isotropic and self- similar in a wave-number range between well separated infra-red m/, m g 
and ultra-violet Mf, M g cut-offs. Specifically 



a i d 2 v e ipx a / mi v \ m f 

f '(., rn„ H^F.}^— ^ « 1 (7) 

with I1 Q/3 the transversal projector in Fourier space, e the Holder exponent and 

G( X , mg ,M g ) = G J^^ Xg (4^) , ^«1 (8) 



2 



, ( 2 7t) 2 p 2ft - 2 A9 \p 2 'M2J ' M 9 
with Holder exponent h. The explicit form of the isotropic cut-off functions Xf an d Xg ls 
unimportant as far as they remain approximately constant for any wave number p in the 
scaling range m,f ,m g <C p <C M/ , M s . It is not restrictive to choose Xf> Xg normalized to 
the unity in the origin. Note that at e, h equal zero the traces of © and ([8]) are proportional 
to the Laplacian of a Dirac 5-function centered at the origin. More generally, the Holder 
exponents e, h determine the spectral composition of the energy injection. For the velocity 
field one finds 

Ie(m f ,M f ) = F a / j^P^Xf , ^ U (9) 

Similarly, the injection for the scalar field is dominated by wave-numbers around the ultra- 
violet (infra-red) cut-off for any h < 2 (h > 2). For the purposes of the present analysis it is 
worth recalling the vorticity representation of Navier-Stokes equation in 2 d: 

(d t + vd x )cu = iyd 2 x uj-- + f ul (10) 

T 

where 

uj = e a pd Xa v p k. f u = e a pd xa f p 

(11) 

ei2 = -e2i = l & en = e 2 2 = 
Equation (fTOl) implies the conservation of the total enstrophy Z 

z = ft, ,^ ( mo,.)^-/a«S,.M^M, (.2) 

whenever the right hand side of (fTOl) is negligible. For power-law forcing, the enstrophy 
injection is also controlled by the Holder exponent e 



M)- 2e 22(0,1), < e < 3 
I z (m f}Mf ) = F I f^P^Xf I ) « < 03) 



d P 4- 2e .. / m / P \ _. J 7 



m 6 r 2£ I Z (1,0), e>3 
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The relations ©, (fT3l show that the energy and enstrophy injections are simultaneously 
concentrated in the ultra-violet and in the infra-red only for e < 2 and e > 3, respectively. 

3. Summary of the renormalization group analysis results 

Let d,A denote the scaling dimension of a physical quantity A. Crudely matching canonical 
dimensions in and © suggests, irrespectively of the spatial dimension, the existence of two 
scaling ranges : a dissipative one for wave-numbers such that non-linear effect are negligible 

d t = 2d x & d v = —d x (l—e) & d e = — d x (l — h) (14) 

and an inertial range corresponding to the requirement of Galilean invariance imposed by 
matching the two terms of the material derivative D t := d t + v ■ d x with the forcing 

d t = d x \2-^-j & d v = -d x (l- & = -/* + !) (15) 

For e , h equal zero the scaling dimensions (fl4l) . (fT5|) coalesce. Physically, the coalescence 
point corresponds to stirring by thermal noise. Coalescence hints at the existence of a marginal 
case in renormalization group sense. This means (see e.g. 1172 1) that scaling dimensions at 
small but finite e , h may be obtained in a Taylor series based at their values in the marginal 
case in analogy to equilibrium critical phenomena where marginality is defined by an upper 
critical dimension specified e.g. by Ginzburg's criterion. A similar scenario seems to apply 
to CQ) for d > 2 and in the absence of large scale friction (r set to infinity). In such a 
case, 11271 l23l l2l fine-tuning the amplitude of the forcing correlation © to be 0(e) yields 
an expansion in powers of e around a Gaussian theory specified by renormalized eddy 
diffusivities vr and kr. Ultra-violet renormalization guarantees that the eddy diffusivities 
are related to the molecular viscosities appearing in (OQ) and © by renormalization constants 

Z v :=— & Z K :=— (16) 

determined at any order in perturbation theory by subtracting all "resonant" terms diverging 
with the ultra-violet cut-offs Mf and M 9 . Technically, this is achieved by identifying the 
ultra-violet divergent part of the one-particle irreducible vertex associated to the response 
functions of the velocity and concentration fields [2J. The result is that, within all order 
accuracy in e, all correlation functions of velocity and concentration fields sampled at well 
separated spatial points scale according to canonical dimensional predictions. In particular, 
the scaling laws (PT41) and (fT5l) correspond to an ultra-violet and an infra-red stable fixed 
point of the renormalization group transformation respectively describing the dissipative and 
inertial ranges. Extending this analysis to the 2d case presents extra difficulties. Already 
in the absence of non-linearities, the correlation function is simultaneously logarithmically 
divergent both in the infra-red and in the ultra-violet. The Ekman term in ([I]) is then needed 
to decouple infra-red degrees of freedom. Once this is done, dimensional analysis shows 
that the renormalization constants (TT6l) are not sufficient to reabsorb all terms divergent 
with the ultra-violet cut-offs in the perturbative solution of (0Q). This is a serious difficulty 
because direct calculations SITU hint that ultra-violet renormalization group transformations 
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using non-local counter-terms may lead to mathematical inconsistencies. In other words, 
renormalization constants should only be associated to coupling constants of local interactions 
in real space. For (0Q) the molecular viscosity is the only coupling constant satisfying such 
requirement [|2l[34l[3]|. The difficulty led to a controversy, summarized in ll34ll . about the very 
possibility to apply renormalization group methods to 2d turbulence. In ll34l it is also argued 
that multiplicative ultra-violet renormalization remains consistent to all orders in perturbation 
theory if the forcing correlation is modified to include a local (analytic) component 

F^(x,m f ,M f ) - F^(x,m f ,M f ) + F^ cal) (x,m f ,M f ) (17) 

with 

From the renormalization group point of view, the replacement is justified by the observation 
that the resulting model generates the same relevant couplings as the original and therefore 
should fall in the same universality class. The merit of (fTTT) is to provide through a a : = 

F (local) 

I v 3 the zeroth order of the extra renormalization constant 

Z a := ^ (19) 

a 

needed to reabsorb all the remaining explicit dependence on Mf in the perturbative expansion 
of correlation functions of the velocity field. In consequence, ll34l predicts that the isotropic 
energy spectrum of the velocity field 

E \M-= f 7^?%-^) f d 2 xe l( *- x ^ v(x,t) ■ v(0,t) y (20) 




(27T) 



admits the expression 



£ M (p) = ^ ('^y'VpMfB Lmt, (»)"•' ) (2d 



v 3 J \ V w 



In (ED the wave number 



V 3 T 3 



Pb OC ^-^r " (22) 



signals whether at small scales dissipation is mainly due to friction (p p b ) or to molecular 
viscosity ( p 3> p b ). The function R has a regular expansion in e for fixed p b . The infra-red 
scaling of (12TT) is then determined by the behavior of R in the limit p J. 0. This limit is inquired 
within the renormalization group formalism by the so-called operator product expansion. 
The outcome ll34l is that the energy spectrum admits the same infra-red asymptotics as in 
dimension higher than two 

E [v] (p) ~ p 1 -^ (23) 

It is worth emphasizing that the above results were derived in If34ll using the vorticity 
representation of the velocity field which holds only in 2 d. As a further check, in 11351133] the 
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same results were recovered by analytic continuation of © in the limit d j 2. In [33 ] the 
analysis extends to the scaling properties of the concentration field and gives 

E [e] (p) ■= J ^SiP-v) J ***** -< 0{x,t)0{0,t) y~ p 1+ ^- 2h (24) 

in agreement with the dimensional prediction (fl~5l) . In summary, according to the 
renormalization group analysis of 041 [351 [331, m 2 d as in 3 d the scaling properties of (OQ), 
([3]) differ for e tending to zero from those of fully developed turbulence. In particular, inverse 
cascade-like scaling is attained only for e equal two and direct cascade-like at e equal three. 
In the following section we will argue that these results are in contradiction with those that a 
phenomenological theory a la Kraichnan would suggest. 

4. Phenomenology a la Kraichnan 

Our starting point are the Karman-Howarth-Monin equation ll30l 

±d^5v»(x,t)\\5v(x,t)\\ 2 y= 

(d t + -< v a {x, t)v a {0, t)y+2u^ d^v a {x, *)0"u a (O, *) y -F\(x) (25) 

and the analogous expression for the scalar field 

U^8v»{x,t) [66(x,t)] 2 y= 

d t ~< [56(x, t)} 2 y+2K^ d„9(x, t)d»6{0, t) y -G(x) (26) 
In (|25T) and (l26l) the notation is 

8v»(x,t) := v> M (x,t)-v fi (0,t) 

60(x, t) := 66(x, t) - 69(0, t) (27) 

The scaling predictions of Kraichnan's theory stem from the asymptotic solution of (|25l) under 
the following three assumptions O: 

(i) velocity correlations are smooth at finite viscosity and exist in the inviscid limit even at 
coinciding points, 

(ii) even in the absence of large scale friction (i.e. r = oo) Galilean invariant functions, and 
in particular structure functions, reach a steady state, 

(Hi) no dissipative anomalies occur for the energy cascade. 

The assumptions (i), (ii) imply that the two point correlation in the absence of Ekman friction 
does not reach a steady state 

-< v a (x,t)v a (0,t) y= Xt- - -< \\5v(x,t)\\ 2 y +. .. (r = oo) (28) 
the constant A being the asymptotic growth-rate. By (Hi) energy dissipation in (|25l) satisfies 

i limlim - limlim I v ~< d^v a (x, t)d^v a (0, t) y= (29) 
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This latter hypothesis is distinctive of two dimensional turbulence: in three and higher 
dimensions the limits are not expected to commute for fully developed turbulence. If the 
bulk of the energy injection Is occurs around a wave-number pf and viscosity and friction are 
such that the adimensional parameter 

K = > 1 (30) 

v 

plays the role of a large Reynolds number then the three hypotheses yield an inverse cascade 
for wave-numbers p in the range p T <C p <C pf with p T = (I £ r 3 ) -1 / 2 and a direct cascade 
for pf <C p <C p T withp T = (z/t) -1 / 2 = p T /lZ l l 2 BSE!- Note that the Kolmogorov scale 
Pk ='■ (-fe/V 3 ) 1 / 4 = p T /TZ 3 ^ is always smaller than the dissipation scale set by the Ekman 
friction. We show below how the same arguments can be adapted to a power law forcing. 



4.1. Inverse cascade: e < 2 

By © and (fT3l) the energy and the enstrophy input are in this case dominated by the ultra- 
violet cut-off Mf. Neglecting nif, the trace of the forcing correlation admits the asymptotic 



expansion (see Appendix A for details) 

{l £ (0,l)- WK^) 2 +,,,} MfX <<; 1 



F a a (x,0,M 



M A f 2 



f) 



jrr(e-l) x4-2 E m f x » 1 



holding for < e < 2. Comparison with the renormalization group results is possible in 
the infra-red region, MfX ^> 1. In order to extricate the corresponding asymptotics of the 
third order structure function, it is convenient to consider first the quasi stationary case for r 
tending to infinity. By hypotheses (i) and (Hi) 10, the asymptotic growth-rate of (1281) in the 
inviscid limit is equal to to the energy injection 

A = Mf- 2£ I £ {0, 1) := M^~ 2£ F* (31) 

By (|25l) the growth-rate sustains the structure function at scales MfX ^> 1 

x wjimn-,*) - nut*-* {i - Jg^gfa + ■ ■ ■} 02) 



In the presence of the Ekman friction (r < oo) (1251) reaches at steady state. In such a case the 
energy injection is balanced by the velocity correlation which far from the infra-red cut-off 

p T = (F M 4_2e r 3 ) -1 / 2 is expected to take the form 

-< v a (x, t)v a (0, t) y= ° 2 f — [l - c x (p T xf 2 + . . . j (33) 

with ci a pure number and (2 to be determined by a self-consistence condition. The 
asymptotics of the structure function acquires a correction 

-< 5v^(x,t)\\5v \\ 2 (x,t) y= 

I (2 + C 2 ) ^ r(e) (Mjx) 4_2e J W 
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Some remarks are in order. The constant flux solution dominates for 

« x « (F M 4 - 2£ r 3 ) 1/2 (35) 

In this range the renormalization group prediction clearly appears as a sub-leading correction. 
Similarly, the two-point correlation adds a further sub-leading term associated to the exponent 
£ 2 - Dimensional considerations yield for £ 2 the value 2/3 whence a —5/3 exponent follows 
for the energy spectrum. The sign of the constant flux term stemming from (1321) . (1341) is 
positive so describing energy transfer to larger scales. The conclusion is of an inverse cascade 
taking place above the forcing ultra-violet cut-off. In such a case, the study of the statistics of 
the passive scalar should recover the results of [[T2l . A priori inspection of (|26l) allows one to 
distinguish at least two sub-cases. 



4.1.1. Scalar field in the inverse cascade and small scale forcing h < 2. For h < 2 
the injection of scalar fluctuations is concentrated in the ultra-violet (thermal stirring). No 
dissipative anomaly is expected. The Karman-Howarth-Monin equation (T26l) yields 

p „/ w™, m2 A 1 - h T(2-h) G x» 

-< *) M*> *)1 ^ nr\h) x ^ (36) 

in the scaling range 

m = {mj^rrig 1 } > x > min {My 1 , M" 1 } := M (37) 

Inferences about the correlation functions of the scalar field can be barely drawn from crude 
dimensional considerations. Accordingly, one has 



-< 



59 Ua 



p 2, 



2 



7-6 h 



E { e](p) — \_i P 3 (38) 

p 3 



which differs from the one stemming from the scaling dimension dg given in (1151) and 
supported by the renormalization group calculations of 11331 . In particular, whilst (fT5l) recovers 
equipartition scaling only for (s,h) = (0,0), (1381) yields a scaling linear in wave-number 
space at h = 2/3. This may indicate the breakdown for h < 2/3 of (1381) and the onset of an 
equipartition-type scaling for the spectrum of the scalar field, see Appendix B for quantitative 
modeling of the phenomenon. 



4.1.2. Scalar field in the inverse cascade and large scale forcing h > 2. In this regime, the 
injection of the scalar field is dominated by the infra-red cut-off. Physically the situation may 
be assimilated to turbulent stirring. The Karman-Howarth-Monin equation for m g x <^ 1 
reduces to 

- G m 4 - 2h {%_ 2 + (m g x) 2h -% + (m 9 x) 2 7h-i + • • •} (39) 



The coefficients of the forcing expansion are specified by the formulae of Appendix A by 



identifying 7 with the function (f> thereby defined. Equation (1391) points at a direct cascade of 
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the scalar field with sub-leading corrections due to the power-law forcing. The dimensional 
prediction for the spectrum of the scalar field is Obukhov-Corssin's 11301 

E [e] (p) ~ p~ 5/3 (40) 
The numerical experiments of |[T2ll support Obukhov-Corssin's scaling in this regime. 



4.2. Local balance: 2 < e < 3 

For 2 < e < 3 the injection of the scalar field is dominated by the infra-red cut-off rrij. The 
quasi- steady state solution for vanishing Ekman friction yields 

A = m 4 - 2£ /£(l,0) :=F o m)- 2£ £ _ 2 (41) 

Correspondingly, the scaling range is set by the condition m^i < 1 alone and (1251) becomes 

^d^5v»(x,t)\\5v(x,t)\\ 2 y= 

\-F mj- 2£ {0 £ _ 2 + (m / x) 2£ - 4 0o + (m / x) 2 e _ 1 + ...} (42) 

As the first two terms on the right hand side of (|42l) cancel out, the third order structure 
function admits the asymptotic expression 

-< 5v»(x,t) \\5v(x,t)\\ 2 y~ -F x**"-* {-^y + (m/X)6 2 2£4 ' 1 + • • j( 43 ) 

with (f) < (see Appendix A). The steady state solution in the presence of Ekman friction 
can be discussed as in subsection 14. II and introduces in such a case only sub-leading terms. 
The cancellation in ([42]> was argued in j6l U\ HI to underlie the direct cascade for standard 
turbulent forcing. The difference here is that power-law forcing dominates the scaling. In this 
range, the solution of (|25T) validates the dimensional prediction (fT31) . The agreement extends 
to the scalar field with two provisos 

(i) h < 2: the threshold for equipartition scaling is h* = e/3 > 2/3: 

(ii) h > 2: the forcing becomes dominated by the infra-red cut-off. Correspondingly, a 
"freezing" of the scaling dimensions at the value for h — 2 may be expected Il2~9l l2~l. 

In summary, the expected spectra are 



p, < h < 



e 



E , (,,) - p 1 -^ & E [e] (p) ~ { y+¥-2^ , I < h < 2 (44) 

-3+- ? 



2s 



p i+ ~ , h > 2 

It should be noted that in this regime for e/3 < h < 2 the predictions of the renormalization 
group and of the phenomenological theory coincide. 



4.3. Direct cascade : e > 3 

For e > 3 both energy © and enstrophy (TT3T) injection are dominated by the infra-red cut-off. 
The analysis of the Karman-Howarth-Monin equation under (z), (ii), (Hi) in the quasi-steady 
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state follows the same lines as in previous subsection 14.21 The second term in the square 
brackets of (|42|) now dominates scaling 



-< 6v"(x, t) \\Sv(x,t)\\ 2 y~ -F x» x 2 m f - 2e { ^ + i^^l 600 + . . . J (45) 

However, (|45T) is just the simplest of the possible scenarios. A detailed analysis of the direct 
cascade stirred by turbulent forcing indicates that scaling in the steady state brought 
about by an Ekman friction may well be characterized by non-universal exponents depending 
upon the value of r. Subsequent numerical investigations ll59l fl4l l68l validate non-universal 
exponents in the steady state. A further inference drawn in Q is the presence of logarithmic 
corrections to the quasi-steady state structure function. For the scopes of the present work it 
is sufficient to observe that the phenomenological theory predicts for e > 3 the "freezing" 
of the scaling dimension d v to a value close to the "naive" direct cascade (d v = d x ). The 
implication for the energy spectra by dimensional arguments is 

{p, < h < 1 

j9 3 - 2/l+ "', 1 < h< 2 (46) 
p- 1+ - , h > 2 

In (|46|) the ". . ." denote non-universal and/or intermittent correction the presence whereof 
is suggested by numerical investigations of the direct cascade stirred by turbulent forcing 



Observation: 

In the above discussion we omitted to discuss the properties of the flow of (QQ), © at scales 
Mf i < 1. The reason for doing so, in the spirit of renormalization group and consistently 
with the numerical experiments of the ensuing section, is that the dynamics is strongly 
suppressed by dissipation effects for scales below Mf. 



5. Numerical experiments 

In order to compare the renormalization group predictions of section [3] with those of the 
phenomenological theory a la Kraichnan expounded in section H we performed numerical 
simulations of the Navier-Stokes equation for the vorticity field (flOl) and the advection- 
diffusion equation for the scalar field (Q~|) with a fully-dealiased pseudo-spectral method [fT9l in 
a doubly periodic square domain of size L = 2n at resolution iV 2 = 1024 2 . Dealiasing cutoff 
is set to k t = N/3. Time evolution was computed by means of a second-order Runge-Kutta 
scheme, with implicit handling of the linear friction and viscous terms. As customary (see 
e.g. B59l[68l) we added to (OQ) an hyperviscous damping (— l) p_1 z/ p _i9 2p f .This is equivalent 
to a Pauli-Villars regularization the use whereof is well justified in renormalization group 
calculations (see e.g. [|72l l34l ). The integration time has been carried out for twenty large 
eddy turn-over times after the velocity fields have reached the stationary state. The stochastic 
forcing is implemented in Fourier space by means of Gaussian, white-in-time noise as in [fT4l 
but with variance determined according to ©. 
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Figure 1. Kinetic energy spectrum for e = 0. Inset: energy flux II^. Parameters values are: 
mf = 1, M f = 341, v-i = 1CT 16 , tvT 1 = 10 2 , F a = 4 • lO" 10 
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Figure 2. Scalar spectra for e = and various values of h. Parameters values are: 
m g = l,M g = 341, k 3 = 1(T 16 , G = 4 • 1CT 10 for ft = 0, m g = l,M g = 341, 
k 3 = 10- 16 , Go = 1.6 • 1(T 5 for h=l,m g = l, M g = 341, n Q = 5 • 10" 4 , G = 2.5 • 1CT 1 
for h = 2.5. Parameters for the velocity field as in Fig.[T] 

In Fig. Q] we show the energy spectrum Er v ] for e equal zero. The numerical spectrum 
exhibits good agreement with the phenomenological theory of section @] with a scaling 
exponent dg = —5/3d p within numerical accuracy. 

Such an exponent is very far from the equipartition-like scaling d £ = d p which is 
the starting point for the renormalization group analysis. The energy flux validates the 
interpretation of the —5/3 spectrum as brought about by an inverse cascade. The energy 
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Figure 3. Kinetic energy spectrum E(k) for e = 2.5. Inset: energy flux He (circles) 
multiplied by minus the unity, and enstrophy flux Hz (triangles). The lines represents the 
injection spectra Ie (dashed line) and Iz (dotted line). Parameters values are: nif = 1, Mf = 
341, v 3 = 10" 17 , r 2 _1 = 10 3 , F = 1 
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Figure 4. Scalar spectra for s = and various values of h. Parameters values are: 

m g = l,M g = 341, k 3 = 1Q- 16 , G = 4 • 1CT 10 for h = 0, m g = l,M g = 341, 
k 3 = 10~ 16 , Go = 1.6 • 1CT 5 for h = l,m g = l, M g = 341, k = 5 • 1CT 4 , G = 2.5 • 1CT 1 
for h = 2.5. Parameters for the velocity field as in Fig. [3] 
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U E {p,m f ) = J ^ J d 2 xe iq ' x -< v a {- q) t){v p d p v a ){x,t) y (47) 

where v a denotes the Fourier transform of v a , is negative and constant in the scaling range 
(see inset of Fig. [B, so signaling the presence of an inverse cascade. 

Breakdown of the marginality assumption at e equal zero is confirmed also by the 




Figure 5. Kinetic energy spectrum for e = 4. Inset: enstrophy flux Tlz- Parameters values 
are: m f = 6, M f = 240, ^ 3 = 10~ 18 , rf 1 = 2, F = 1 



concentration spectra i£[0] shown in Fig. |2l at h = 1 a spectrum E[q] ~ p 1 / 3 well fits the 
numerical results, which definitely rule out the E\g\ ~ p^ 1 prediction of the renormalization 
group theory (|24)) . Furthermore and in agreement with sub-section 14. II the spectrum of the 
scalar field undergoes a transition at h = 2. Above that value the scaling exponent freeze 
to the value —5/3 corresponding to a direct cascade, in agreement with the results of [fT2ll . 
Fig. |2] illustrates the phenomenon for h equal 2.5. It should be noted that for h equal zero 
an equipartition-type scaling (i.e. linear in wave-number) is observed for the spectrum of the 
scalar field. It is worth stressing that this is not sufficient to infer equipartition of the full 
statistics of the scalar field. It is known in analytically tractable cases of turbulent advection 
of scalar fields 11281 I2T1 l22l |54| that equipartition-type scaling of the two-point correlation may 
well co-exist with highly intermittent statistics even in the decay range of a scalar field. 

In agreement with the phenomenological theory, for 2 < e < 3 no cascade is 
observed. Fig. [3] illustrates the situation at e = 2.5. The energy injection spectrum, defined 
as = dpp 3 ~ 2e Xf y-pf , jjij is dominated by IR contributions, while the enstrophy 

injection spectrum Iz{k) = j Q dpp 5 ~ 2e Xf f —f , j is still ultraviolet divergent (see inset of 
Fig. [3]). In this situation the steady state is characterized by a scale-by-scale balance between 
the fluxes and the injection spectra. The resulting energy spectrum scales as E^{p) ~ p 1 ^^ 
(see Fig. [3]). The spectra for the passive scalar shown in Fig. |4] are in agreement with the 
prediction (l44j) . 

Finally, Fig J5] evinces the onset of a direct cascade for e > 3. The enstrophy flux, 

IL z (k 1 m f )= [ 7^%9f? [d 2 xe iq - x ^ u(-q,t)(v%u)(x,t) y (48) 

Jm f \£ 7r ) J 

is approximately constant and positive (indicating transfer towards smaller spatial scales) in 
the numerically resolved scaling range (see inset of FigJ5]). It should be emphasized that 
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numerical evidences of ||59l [681 [T4l uphold non-universal dependence of the kinetic energy 
spectrum upon the Ekman friction in the direct cascade regime. Since the analysis of such 
effects lays beyond the scopes of the present work we replaced for e > 2 the Ekman 
friction with an hypo-dissipative term (— l) q+1 T^ 1 d^ 2q v which is expected to suppress the 
aforementioned non-universal corrections to the spectrum [48] . 

6. Discussions and Conclusions 

Our numerical experiments support without possible ambiguity the scenario set by 
Kraichnan's phenomenological theory. The physically relevant order parameters to describe 
the qualitative behavior of (OQ), © are the total energy © and enstrophy (TT3T) injections. 
Whenever they coherently act on small (large) scales an inverse (direct) cascade is observed. 
In the intermediate case 2 < e < 3 no cascade takes place and a local balance scaling 
takes place. The renormalization group scaling exponents are for e , h < 2 at most a sub- 
leading correction to the inverse cascade scaling. It remains to be clarified the origin of the 
discrepancy between the two theories. The renormalization group analysis of Oil [33l E21 is 
very thorough and satisfies all the self-consistency requirements that in the theory of critical 
phenomena are known to produce correct scaling predictions. Thus any trivial explanation 
of the discrepancy can be safely ruled out. The indication of the existence of a constant- 
flux solution stems from the Karman-Howarth-Monin equation. In the renormalization group 
language the relations closest to (|25T) are those satisfied by the composite operator algebra 
which includes the energy dissipation operator (see e.g. section 2.2 of [2]). One may 
speculate that non-locality of forcing plays a role different than in 3d in determining the 
scaling dimensions of the elements of the operator algebra. It is however difficult to see how to 
consistently formalize this observation. One scenario that deserves to be further investigated 
is in our opinion the following. Perturbative renormalization implies the assumption that 
the infra-red stable fixed point governing the scaling regime emerges from a bifurcation at 
marginality from a Gaussian fixed point. Such an assumption is usually verified in critical 
phenomena but is not a necessary consequence of a general non-perturbative theory. In 
particular, there are examples of field theories where it is possible to give evidence that 
scaling is dominated by a fixed-point emerging at marginality from bifurcations from non- 
perturbative, non-Gaussian fixed points. A concrete case is discussed in ItBTl 1501 Q. Here 
a model of wetting transition indicates a scenario which could apply also to 2<i-turbulence. 
The scaling predictions associated to an infra-red stable fixed point captured by perturbative 
renormalization group analysis is numerically seen to be dominated by those associated to 
a second fixed point not bifurcating from the Gaussian fixed point at marginality. The 
existence of this fixed point can be exhibited only by a non-perturbative construction of 
the renormalization group transformation. The price to pay is however the introduction of 
truncations of the Wilson recursion scheme which cannot be a-priori fully justified. If we 
accept such point of view, the existence of the Karman-Howarth-Monin equation should be 
interpreted as an a-priori indication of the existence of a non-perturbative fixed point. 

X We thank L. Peliti for drawing our attention to this point and to refs. lf5Tl l50l . 
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Appendix A. Asymptotics of the forcing correlation 

Consider the scalar correlation 

with 

X(0,0) = 1 (A.2) 

In order to extricate the asymptotics in the range m <^ p <^ M two cases should be 
distinguished depending upon the sign of r\. 

• If r\ < the integral is infra-red and ultra-violet convergent at finite point separations in 
the absence of cut-offs: 

F(x,m,M) ~ F / , ., = — ^-r K2; 21 (A.3) 

with Vl d the solid angle in <i-dimensions. 

• If T] > dA. II) is convergent only if the infra-red cut-off is retained: 

F(x,m,M)c±F(x,m,oo) = F J ^_ X ^_, J (A.4) 

The integral can be estimated by inverting its Mellin representation with help of Cauchy 
theorem (see e.g. 11561 ) 

F(x,m,w) = F ^^^( mx )2C0( C); + ^ 

JStC-%00 {^m) 2 

The holomorphic function 



(2 Tr)^ 2 1+2 C+^ r ( d+2 2 +V ) Jo w w " 

is by hypothesis analytic at least in a stripe for < —rj/2. Furthermore, by (IA.2I) . 
the integral in (IA.6I ) generates a simple pole for £ equal zero. Cauchy theorem yields for 
m x <C 1 the asymptotics 

F(x, m, oo) = F m~ r? |<^_2 + (ms)' ? </> + (mx) 2 0!_h + . . • j (A.7) 

having used the notation 

a :=-lim(C-a)0(C) (A.8) 
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Appendix B. Large scale zero-mode and power-law forcing in the Kraichnan model of 
advection of a concentration field 

We refer the readers for definition and details on the large scale decay properties of the 
Kraichnan model to ll28ll2Tll22ll54l . The energy spectrum of the Kraichnan model is exactly 
known. In 2 d it takes, modulo irrelevant constant factors, for isotropic forcing the form 

dp pBi{p p) r p da 
p k + D p& J a 



E[e]{p) = I I -<r 2 G(*,m g ,M g ) (B.l) 



with B v (x) the Bessel function of order t], D the eddy diffusivity of the advecting velocity 
field and G given by ([8]). The Holder exponent £ is a free parameter in the model. Turbulent 
advection corresponds to £ equal 4/3. From (IB.ll) it is straightforward to check that 

roo dp p 2 Go _ (K\^ Gp-Kp 1 I 

"Jo p K+Dpip2- 2h - \D) K £ sin(i|2) ' 2 

E[ 9 ](p) ~ { (B.2) 

n l+£-2A f°° dp Go Bi(p) ____ GopW 2fe r(fe-§) 7 i 

w Jo p zv+e-a* - 2 2+c-2"Dr(2-fe+|) ' ^ 2 

G is a dimensional constant, the value of which is irrelevant for the present considerations. 
Setting £ = 4/3 the results of subsection (14.11 ) of the main text are recovered. The direct 
cascade results (l46l) are recovered instead by setting £ = 2. The asymptotics (IB.2I) hold 
under the assumption of infinite integral scale 1 of of the spatial forcing correlation. In the 
language of I12T1I221 this means that the "charge" G(0, 0, M g ) is vanishing. For h > 2 the 
forcing becomes infra-red dominated. At scales m g x 1 



E [9] (p)~G(0,m g ,oo)p^ f*P(M?M. (B.3) 

Jo P u 

whilst in the opposite range m g x 3> 1, for (5(0, m g , oo) > 



E m (p) ~ G(0, mff) oo)^ / ^ p 1 ^^ (B.4) 
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